Kondo and charge fluctuation resistivity due to Anderson impurities in graphene 
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Motivated by experiments on ion irradiated graphene, we compute the resistivity of graphene 
with dilute impurities. In the local moment regime we employ the perturbation theory up to third 
order in the exchange coupling to determine the behavior at high temperatures within the Kondo 
model. Resistivity due to charge fluctuations is obtained within the mean field approach on the 
Anderson impurity model. Due to the linear spectrum of the graphene the Kondo behavior is shown 
to depend on the gate voltage applied. The location of the impurity on the graphene sheet is an 
important variable determining its effect on the Kondo scale and resitivity. Our results show that for 
chemical potential near the node the charge fluctuations is responsible for the observed temperature 
dependence of resistivity while away from the node the spin fluctuations take over. Quantitative 
agreement with experimental data is achieved if the energy of the impurity level varies linearly with 
the chemical potential. 

PACS numbers: 



I. INTRODUCTION 

A logarithmic upturn in the resistivity at low temper- 
ature has been observed in graphene with vacancies^. A 
fit to the temperature dependence of resistivity with con- 
ventional Kondo effect yields a large Kondo temperature 
(with Tk ~ 30 ~ 90K) which shows a non-monotonic 
behavior with respect to the gate voltage^. The vacan- 
cies in the graphene sheets are induced by ion irradiation 
in ultra-high vacuum and the magnetism in sputtered 
graphite has been experimentally observed 2 - - — . Our goal 
is to study whether Kondo effect alone in graphene can 
explain the experimental results in Refill. 

We start with the Anderson impurity model^ to study 
the impurities effect on transport. In the local moment 
regime where the impurity occupation for a given spin 
Ud :S — 0.5 we use Schrieffer Wolff transformation^ to 
write down the Kondo model from the Anderson impu- 
rity Hamiltonian. Since we are interested in the resistiv- 
ity due to impurity spin fluctuations we study the Kondo 
model by standard perturbation method. The perturba- 
tive approach fails at T ~ Tk with Tk representing Kondo 
temperature but works for Tk <C T. We compute the 
scattering rates in this weak coupling regime. The scat- 
tering rates are determined via perturbative calculations 
of the T-matrb*i2r— . Kondo effect in the pseudogap sys- 
tem has been explored in the context of graphene as well 
as in that of d-wave superconductor— ~— via various dif- 
ferent approach such as NRG or mean field approach^. 
The advantage of our approach is the ability to determine 
the high temperature behavior of the scattering rate and 
resistivity accurately within perturbation theory. 

We assume a dilute concentration of impurities and 
ignore the spin spin interactions such as the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) interaction. In graphene 
these interactions, in addition to being oscillatory with 
distance between impurities, depend on the sublattice on 
which the impurities are locate d 21 ' 22 . For chemical po- 
tential at the Dirac point our results are in agreement 
with the prediction of the existence of an intermediate 



coupling fixed pointi£~— . Near the node the exchange 
coupling J needs to be larger than a critical value J c 
to have the Kondo effect. The dependence of Tk on the 
chemical potential is qualitatively different for /i <C Tk 
and fj, 3> Tk- For impurities breaking the lattice sym- 
metry, a power law in T divergence of scattering rate is 
obtained for /i <C Jfc while a logarithmic divergence ap- 
pears for > T|j. For impurities preserving the lattice 
symmetry, a power law in T 3 divergence of scattering rate 
is obtained for /i -C while a similar logarithmic diver- 
gence appears for /i 3> Tk- For both cases the scaling of 
resistivity with single Kondo temperature breaks down 
in the vicinity of the Dirac point. Our results for Kondo 
temperature obtained within the T-matrix formalism is 
in agreement with the mean field results for the devel- 
opment of the Kondo phase^iZ. The resistivity obtained 
displays different chemical potential dependences. For 
impurities breaking the lattice symmetry the resistivity 
decreases as chemical potential increases while for impu- 
rities preserving the lattice symmetry the resistivity in- 
creases as chemical potential increases. For the same set 
of physical parameters the dominant source of resistivity 
is from impurities which breaks the lattice symmetry. 

We also explore the region near the empty orbital to 
mixed valence one in the Anderson impurity model to 
find the resistivity due to charge fluctuations. From the 
numerical RG 14 -^ the Kondo effect is suppressed as the 
critical exchange coupling J c — > oo for chemical potential 
close to the Dirac point. Thus we use unrestricted Hatree 
Fock method 2 ^ on the Anderson impurity model to find 
the resistivity near the empty orbital regime. The re- 
sulting resistivity shows similar dependence on chemical 
potential as well as dominance from symmetry break- 
ing impurities as the resistivity obtained in the Kondo 
model. Near the node the Kondo scale, extracted from 
the logarithmic temperature dependence region on re- 
sistivity, yields a Kondo temperature comparable to the 
observations in the experiment in Refill while away from 
the Dirac point the extracted Kondo scale is higher than 
experiment by one order of magnitude. 
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By combining the charge fluctuation effect for [i ~ 
and Kondo effect (spin fiuctuations) for finite [i we ob- 
tain Kondo temperature dependence on fi qualitatively 
consistent with experimental results 1 with gate voltage 
less than 301^. Our conclusion is that the observed ex- 
perimental results, albeit fitted well by Numerical RG 
for conventional metal Kondo model 24 , cannot be solely 
explained by Kondo screening in all range of chemical po- 
tential. For chemical potential near the node the charge 
fluctuations is responsible for the observed resistivity 
temperature dependence while away from the node the 
spin fluctuations take over. 

This article is organized as following: We start with 
the Anderson impurity Hamiltonian to describe dilute 
impurities physics in the graphene system. To study the 
local moment regime we use Schrieffer Wolff transforma- 
tion to obtain Kondo model from Anderson Hamiltonian. 
In section 3 we evaluate resistivity due to spin fluctua- 
tions, with different impurity locations, by perturbation 
computations on the Kondo model. In section 4 we com- 
pute resistivity due to charge fluctuations when impurity 
occupation is close to zero by using mean field approach 
on the Anderson model. In section 5 we show numeri- 
cal results of temperature dependence of the resistivity 
with different symmetry and mechanism. In section 6 we 
compare our results with the experiment in Ref[l|. The 
results are summarized in section 7. Two appendixes con- 
tain derivations for the perturbative results in the Kondo 
model. 

II. HAMILTONIAN 

We start from graphene Hamiltonian in the presence 
of dilute impurities described by the Anderson impurity 
Hamiltonian^ 

H = Hg + H d + Hu + H v (1) 
H g = -t^2 faa\ s b k , s + h.c. - n(al s a k , s + &fc, s &fc,s) 

k . S 

H d = ^ e dd\d s 

S 

H v = uS s d 8 S_ s d- s 

H v = X>«%4,. + VbAjds + h.c. 

H g is the nearest hopping in the momentum space with 
t ~ 2.7 'eV being the nearest neighbor hopping strength. 
fi defines the Fermi level measured from the Dirac point, 
at and tit are the particle creation operators on the 

a and b sublattices. fa — ^2 i=1 e %k ' ai with o*i = a a x, 
&i = o-q{^x/2 + v3j7/2), and CI3 = ao(—x/2 — y3y/2) 
being the nearest neighbor lattice vector, ao ~ 1.42 A 
is the lattice constant. Hy describes the hybridization 
between the impurity level and graphene electrons with 
Vk,a/b = Si=i Vi,a/be tk ' ai ■ Hjj describes the Coulomb re- 



pulsion on the impurity level and Hd is the Hamiltonian 
describing the localized level of d electron. We diagonal- 
ize H g by defining c\ s± = {a\ s ± (fa/\fa\)ti[ s ). In this 
basis the H g term becomes 

H 9 = ^2int\fa\ - n)c{ sn c ks n (2) 

ksn 

with n = ± denoting the conduction and valence bands. 
The hybridization term Hy in this rotated basis is 

n=± k,s 

with Q kn = (V k>b +nV* a <t>* k /\fa\)/(V2V). Denote e k>n = 
nt\fa\ — fi as the energy of the bands evaluated from 
chemical potential /j, and {k} = (ksn) as combinations 
of momentum, spin, and n the band index. The Ander- 
son impurity Hamiltonian describing the impurity in the 
graphene can be written as 

H = ekncl sn c ksn + ^2 £dd\d s + U dld s dL 8 d- s 
{k} 

+ VYpkncl sn d s + QlJ s Cksn] (3) 
{fe} 

To explore the local moment regime where impurity oc- 
cupation for a given spin rid,s — 0.5 we perform Schrieffer 
Wolf transformation to project out the charge degree of 
freedom!. The exchange Hamiltonian or Kondo model 
obtained after this transformation with the additional 
term Hi nt oc S r ■ S r i describing spin spin interaction at 
different sites is given by H = Hg + Hi mp + Hi nt with 

Hg = £knC ksn Cksn (4) 

{fe} 

Himp — QknQk'm(K5 s >s — ■ (Js's)c k , s , m Cksn 

{fc},{fc'} 

Him = -^W{r-r')Sr-S r - 

r,r' 

Here J ~ V 2 (l/(e d - n) - l/(e d + U - fi)) and K ~ 
{V 2 /2){l/(u - e d ) + l/{u - e d - U)). The interaction 
between impurity spins H int are added for the inclusion 
of spin spin interaction but is assumed to be small due 
to small concentration of the impurities in this article. 
Including this term would lead to a time dependent im- 
purity spin via S(t) = e H, - ntT Se~ HintT with r being the 
imaginary timei^. 

For impurities preserving the C% v point group symme- 
try of the triangular sublattice in the graphene system 
the factor \Qkn \ oc \fa\ while for impurities breaking the 
symmetry the factor |0fc n | is a constant. To evaluate the 
resistivity due to spin fluctuations we use Eq.Q as the 
starting Hamiltonian. We use perturbation expansion on 
the one particle Green function's T-matrix to compute 
scattering rate and from Boltzmann transport to obtain 
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linear response resistivity in both impurity breaking and 
preserving the lattice symmetry cases. We use mean 
field approach on the Anderson impurity model shown 
in Eq.([3]) to obtain resistivity due to charge fluctuations 
in both symmetry breaking and preserving case. The fol- 
lowing two sections are the computation results for each 
cases mentioned above. 



III. RESISTIVITY DUE TO IMPURITY SPIN 
FLUCTUATIONS 

To study the resistivity due to spin fluctuations we 
start with the Kondo Hamiltonian shown in Eq.(0}. We 
calculate transport properties from the T-matrix which 
is related to single particle Green's function byi£ 

G{k'},{k},a', a (iui,iu2) = G{ k >},{k}( iuj i) (5) 
+ G \k'}Ak'}^) T W}Ak}-a'A i ^i^^)G\ k}{k} {iu 2 ) 



where G° fc , } {fc} (zu;i) 



5k',kSs>, s 5m,n{iwi - £fe,n) 1 and 



a, a' are impurity spin state indices. This expression is 
related to time dependent Green's function by 



G{k>},{k},u' ,u{iui, iuj 2 ) 
rP rP 



(G) 



o 



with id — (2r + 1)ttT = (2r + 1)tt//3 and r being inte- 
gers (We put Boltzmann constant fee = 1 to simplify 
the notation). Based on perturbation in Hi nt this time 
dependent Green's function can be written as 



G m , {kha >ATy) = ^zr 1 J2-n\ dn --- dTn E E E ( 7 ) 

n=0 Qfl)gi)",9m,9n si,s' x ,..,s n ,s' n ni,irii,..,n m ,m n 

[ T s(Q*k 1 .n 1 Qk' 1 , mi { KS s' 1 ,s 1 ~ ^( T l) ' Vs{, Sl ) ■ ■ ■ Q k n ,n n Q k' n ,m n {K5 s ' ntS „ - ^S(r n ) ■ CT a ;, Sn ))]a',a 
(T T (cfe' s ' m (r)cfe sn (T )Cfe^ s ^ mi (ri)cfe 1 s 1 n 1 (n) . . . Ck' n s' n m n (T n )ck n s n n n (T n ))H e 



Here T$ and T T are the time ordering operators and (Z) 
is the S-matrix. The first order in J and K is given by 

- ^(T s >sS a ' a {i{^\ - CJ 2 ))] (8) 

Here S(iw') = Jjf dre iu ' T S{T) with J = 2nrk B T. For 
Hrkky = we can simplify above expression by noting 

I 



that S{ud') — /3S6u'.o and we get 

T '{n{fe}, a ', a ( ia 'i' iw 2) = ^ 1 , U3 r {i)}, { fe}, c< ', a K) ( 9 ) 

withT {k'},{k}, a >, a ( iuj i) =®*kn e k'm(K-i°S>sS a > a ). The 
general second order result of T-matrix, with n = 
l/(e**.« + 1) and = iV" 1 E fc , n |e fc „| 2 (z - e^)" 1 , 
is expressed as 



To focus on the Kondo contribution to the scattering 
rate we may set K — and we set Hi„ t = by 
assuming dilute impurities. For RKKY type of spin 
spin interactions the interaction strength decays as 1/i? 3 
for symmetry breaking or 1/i? 7 for symmetry preserv- 



ing case 7 -. Thus for sufficient dilute impurities we may 
treat Hi nt = 0. In this limit Eq. tjTU]) is simplified to 
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T, 



fk>},{k}, a >, a ( lUJ i,iu 2 ) = ^ Wl ,o, 2 r { ^ }i{fe}iCt , !a (iwi) wi 



(2) 



with 



f^W = ©L0 fe ' m (^) 2 [5(5 + 1)F(*)£ S ,, S 5, 



and VK- 2 ) is the digamma function. Analytic forms can 
be obtained in two asymptotic limits by using the asymp- 
totic forms of the digamma function. For /3fi <C 2tt we 
have 



tanh( — —)(S a 'a ■ o- s i s )] 



(11) 



q,l 



For noninteracting spins the third order perturbation, af- 
ter taking a trace over conduction electron spins and ap- 
proximating the reduction of three-spin correlation func- 
tions to two-spin correlation functions 10 , is given by 



y, |e g2 „ 2 | 2 tanh(^^) 



9i,ni 



|e 9ini | 2 



92, "2 



c «3l.ni 



c <?2,™2 



(12) 



From Eq. (fTTj) and Eq. (fT2|) we may define a general func- 
tion R(z) which we need to evaluate in computing the 
T-matrix 



q,n 4 ' 



,anh(^) 



(13) 



By using e q , n = nt\(f) q \ - /j, ~ n^&\q\ - fi = n\e\ - fi we 
may write the continuous form of Eq. (|13[) as 



R(z) 



9nt 2 



A ie i 2 

A e - (z + MJ 



tanh 



(«^)(») 



where A is the linear spectrum cutoff. In the following we 
will separate the discussions into two cases^: One with 
impurity interactions breaking the C% v lattice symmetry, 
in which |0fc n | 2 = |, and we denote R(z) — Rsb( z ) m 
this case. Another with impurity interactions preserving 

the symmetry, in which |6fc„| 2 = 9 ^ fc g a ° 
denote R(z) = Rsp(z) in this case. 



^jt, and we 



A. Cs v symmetry breaking impurities 

For the case of symmetry breaking the cutoff scheme 
we choose for a linear density of states with a cutoff A is 
multiplying the argument of right hand side of Eq.(fl"4|) by 
t 2^_ A 2 and extend the integration limit from ±A to ±oo. 
The resulting Rsb{z), with details shown in Appendix. 

El is 



Rsb(z) 



-I 



Ar 



9?rf 2 I 7T 

xdx 

-F{x,(J,,; 



P 



xdx 
x - A 



F(x,fj,,z) 



(15) 



jo X + A 
Here F{x, fi, z) is defined as 
Ml + (M + 



F(x,fi,z) 



x + + z) 



.0z 2A 2 (z + /x) ■ 
*2^(z + ^) 2 + A 2 . 



a; — i(/i + z) 



Rsb{z 



ttA 



4 7 + 41n(2) -4-3C(2) 



9?rf 2 

(4 7 -2C(2) + 8)(z + /i) 



/3 



/3z(z + n) 



ln(l + ( 



/8(« + /i) 



tan _1 ( 



) 2 ) 



p{z + fx) 



) 



2A 2 (z + M ) I .$z s 
> + /z) 2 +A 2Vl 2 l 2n' 



(16) 



where 7 ~ 0.577 is the Euler constant and f(2) = 7r 2 /6 
is the Riemann zeta function evaluated at 2. In the limit 
2ir we get 



Rsb(z) 



9irt 2 



)(z + M ) tan-^^M 



2tt 



2 + /x 



-7rA - 

-|z + /i|ln(l 



I Ml 



SH^-tan"^ 

7T 2 Z + /i 



)-l) 



/i 2 + z 2 ' 



2A 2 (z + ,u) 1 _ 
(z + At ) 2 +A 2 ^ ( 2 Z 27r y 



(17) 



In Eq. tTTrJl) and Eq. tfTT]) we have assumed < z < 
fx. Using the Boltzmann equation with relaxation time 
hypothesis 2 ^ and noticing that the honeycomb symmetry 
is broken by the impurity we find the scattering rate is 
related to the T-matrix by 



Tss(efe,ra) 



TTTlj 

z ~h 

nni 



S(e k:n - efc/, m )|Tfc,fc'| 2 (l - cos6k,k') 
S(ek, n - efe',m)|2fc,fc'| 2 (l - cos6» feife * 



dk' 
die' 



= 7—-3[J-k,k(£k,n)\ 

n 



(27T) 2 

(18) 



The second line of Eq. (fT8"|) represents the scattering pro- 
cess related to different Dirac points in the Brillouin zone 
and the third line is the scattering event within the same 
Dirac cone. We have used the fact that Q[T k k (e + )} = 
-7r£ fe ' \T k , k >(e + )\ 2 5(e-e k ,) and \T k , k ,(e+)\ 2 is'indepen- 
dent of angle 9 k ^ k i between momenta k and k' in the sym- 
metry breaking case in the above equation. The scatter- 
ing rate, Tq B {u) with tsb(w) being the relaxation time, 
to third order is 

tvrg&u) = 4n / 5(5 + l)[(A)(|) 2 |c + M | 



\u + (19) 
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T(eV) //[eV] 



FIG. 1: The function Rsb(0) plotted as a function of temper- 
ature fcsT in the unit of eV. Energy cutoff A = 3.5eV and 
t = 2.7eV. 



FIG. 2: The Kondo temperature as a function of chemical 
potential for various values of J. Energy cutoff A = 3.5eV 
and t = 2.7eV gives | J c \ ~ 4.68eF. 



The Kondo effect is reflected in the divergence of the 
relaxation time in the parquet approximation. This in- 
volves treating the cubic term as the first in an infinite 
series which is summed to give 



^» = 4n z 5(5 + l)(^)(-) T 



/'I 



m[R(uj)} 



(20) 



We are mainly interested in the DC response so we shall 
study the relaxation time when uj — > 0. FigJT] shows the 
function i?ss(0) plotted as a function of temperature for 
different chemical potential. We define the Kondo tem- 
perature as the temperature when the relaxation time 
diverges when lu — > 0. For the case <C 2tt the singu- 
larities from 1 + J%t[R(0)] ~ can be expressed, by using 
Eq.QI]), as 



4J 



97rf 2 



7rA 



4 7 + 41n(2)-4-3C(2) 



fa 

7rA 



4 7 + 41n(2) -4-3C(2) 



-1 

J 



(21) 



where J c = — 9i 2 /4A and fa = l/T k . Thus for chemical 
potential /i « T we have no Kondo effect if |J| < |J C |. 
As one increases the chemical potential we may include 
the linear order of /j, in Eq. (|16|) and obtain the expression 
for Kondo temperature as 



n = - 



7rA(l-^) + (41n(2) + C(2)-4)/i 
4 7 + 41n(2) - 4-3C(2) 



(22) 



f ) f and V(^) 



where we have used tan 1 {-r l —. — , ., ...... . , 

— 7— 2 ln(2). Thus the Kondo temperature increases with 
increasing chemical potential. 

In the opposite limit where (3[i 3> 27r we use Eq. fTT)) 
to obtain the Kondo temperature as 



T k = ci/xexp 



(23) 



where ci = exp( 7 + ln(2) - 1 - 4/7r)/27r ~ 0.058. Eq. 
can be expressed as T k cx exp(( J— J c )/(p g J) with p g cx fi 
is the electron density of states of graphene. Com- 
pared with the Kondo temperature in conventional metal 
Tfc cx exp(l/p g J) there exists critical value of exchange 
coupling for Kondo effect to be realized in this two dimen- 
sional pseudo gap system. Fig[2] shows the Kondo tem- 
perature as a function of the chemical potential for var- 
ious values of J. For J/J c smaller than 0.87, the Kondo 
temperature is smaller than the chemical potential for the 
range shown. In this regime Tk is given by Eq. (|23[) and 
is exponentially smaller than the energy scale set by the 
chemical potential. As J / J c approaches one from below, 
the Kondo temperature grows faster than the chemical 
potential. As Tk(/J.) ~ /i/27r the exponential behavior 
cross over to the linear dependence shown in Eq. ([2"2l . 

Given the relaxation time we obtain the linear response 
conductivity as 

2e 2 ^ f o . - df dk 



de k , n (2tt) 2 (24) 



2 A„,2 (-00 



2e 2 4u 



3 9na 



2 roc 

fla / de|e|r SB (e-/x) 

l J -00 



tk 



For small chemical potential or f3fi <C 2n we use 
Eq. (ll"6|) and Eq. ([22|) and approximate / g(e)(— df/de) ~ 
g(e)/(2T). We get the resistivity at low chemical 
potential as 

„ . . 3 nS(S + l)alm ( J\ 2 9nt 2 ( 1 
PsbK 1 ) - 7^2 



2e 2 Hv 2 F \2 J 2J \r(T,fj,) 

SimrSiS + l) fJ\ h ( 1 



with 



2 J e 2 \r(T,fi) 



-(T,/i) = (4 7 + 41n(2)- 4- 3C(2)) (T-T fe ) 
- (2 7 -C(2)+4-2^(l))^^ 



(25) 



6 



For temperature fj,/2n < T < /i but higher than T k the 
same approximation scheme gives 



3nn I S{S + 1) 



r(T, /i) = (47 + 4 ln(2) - 4 - 3C(2))(T - T fc ) (26) 

Thus we see that for T > /j, the Kondo contribution to 
resistance is not determined by a single scale T k . For 
temperature ranged between [ijlix < T < fi the scaling 
of the resistivity goes like 1/(T — T/.). This power law 
behavior indicates that at sufficient low chemical poten- 
tial the magnetic impurities are not completely quenched 
while a logarithmic behavior is expected in the conven- 
tional metal case. 

For large fi or j3fi 3> 2tt a Kondo effect similar to 
magnetic impurities in the conventional metals is ob- 
tained. For large chemical potential we approximate 
df/de ~ — 8{e — fi). Under this approximation the re- 
sistivity p s SB {T) = a s^ T) is given by 



f 



Psb(T) 



3 irS(S+l)alm 



2e 2 



hv 2 F 



i + m(R SB (o)) 



(27) 



Use Eq.([Tf|) for K(i? SB (0)) with T > T k we get 



2e 2 



2 J 2J/z 



4 \2jl) e 2 



T 



T 
i 



(28) 



B. Cs v symmetry preserving impurities 

For the case of impurities preserving the symmetry of 
honeycomb lattice the cutoff scheme we choose for a lin- 
ear density of states with a cutoff A is multiplying the 
argument of right hand side of Ea.(fI3|l by A 4 /(e 4 + A 4 ) 
and extend the integration limit from ±A to ±oo. The 
resulting Rsp(z), with details shown in Appendix iBl is 



4 r2A 4 f °° 

Rsp(z) = i 



dx- 



9?rf 4 I 7T J x 4 + A 4 

(z + /i) 3 A 4 1 .f3z 



F(x,fi,z) 



/3/j,/2ir < 1 we have 



Rsp(z) 



1 - 



f tt (7r-41n(^))^ 



~ 9^t 4 I V2 
; 4 + 81n(f))z 



7rA 

/ 4(l + 21n(2)) 
9tt 



^ + ^ 1 




2z/i 
lY 2 " 



7T 

/3A 



+o((^) 2 ,(|) 2 )} 

For the opposite limit (3^l/2tt > 1 we have 



(30) 



Rsp(z) 



4A 3 



f 7T 2tt {z + h) (2h /o fi 2 

97ri 4 lV2 /3A + ttA I A 1 + j A 2 



+ 2( 1 + ^)tan^)) + -^ (7r - 
+^3~ ( - M( z + M) 2 + ( z + A 1 ) 3 tan ~ 



4h(i)) 



2 + \l 



+£) In f it 
3 / V 2tt 



m 2 (m 



(/j + ^) 3 
A 3 



In 



A 3 

(a^ + ^) 2 

/z 2 + (M + *) : 



0} 



(31) 



Similar to Eq.(Tr6"|) and Eq. ([T7|) we have assumed < 
z < fi. Using Eg. (I18p but with the appropriate re- 
laxation times determined in this section, we com- 
pute the resistance. For this case ^s[T k k (e + )] = 
-7r£ fe ' \T k ,k,(e+)\ 2 5(e-e k ') and |T fe , fe ,(e+)| 2 is indepen- 
dent of angle 9 k . k > between momenta k and k' since 
\T k , k >(e + )\ 2 cx |/c| 2 |/c'| 2 . The scattering rate to third order 
is 



f\T, 



iP (u) = 4 ni S(S + 1) 



,J_ V J 



\$t 4 )\2) ' /if 
\lu + ^[Rspicu)]] (32) 



The expression for the relaxation time t$p, within the 
same approach as the previous section, is 



hTgp(u) = 4njS(S+l) 



1 + J5R[i? S pH] 



03) 



A* + (z + ^2-'y (29) 
Ae l3 f - (z + n) Ae 1 ^ - (z + n) 



The DC conductivity is related to the relaxation time 
with uj — > 0. Fig[3] shows the function Rsp(0) plotted as 
a function of temperature for different chemical poten- 
tial. Rs p (0) shows small variations with temperature ex- 
cept when temperature is close to zero where exponential 
growth with decreasing temperature is observed. For the 
case j3fj, <C 2tt the singularities from 1 + Jdi[Rsp(0)] ~ 
can be expressed , by using Eq. (|30| . as 




Analytic forms of Rgp(z) is obtained by taking the 
asymptotic behavior of digamma function in the fol- 
lowing two limits: Ph/2it <C 1 and /3fi/2n 3> 1. For 
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FIG. 3: The function Rsp(0) plotted as a function of tem- 
perature fcflT in the unit of eV. Energy cutoff A = 3.5eV 
and t = 2.7eV. Compared with Fig[2] Rsp(0) shows more 
temperature variations when T — > 0, indicating the Kondo 
effect can only be observed on lower temperature compared 
with symmetry breaking case. 
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Fig|4] shows the Kondo temperature as a function of 
chemical potential for various exchange coupling strength 
J. For J/J c smaller than 0.75 the Kondo temperature is 
always smaller than the chemical potential for the range 
shown. J/J c approaches one from below and for T k — 
fi/27: the exponential dependence on /i crosses over to a 
power law. 

Given the relaxation time we obtain the linear response 
conductivity as 



'SP 



(T) = 



2e2 
3 

2e 2 



n 

Av 2 



df dk 



3 97ra 



de Kn (2tt) 2 

' r aw < u - f 



(36) 



For ^ < 2tt, T > (i, and T > T k we use Eq.(|5D|> 
and Eq. (p4|) and again approximate f g(e)(— df/de) ~ 

!a-T 3( £ )/( 2T )- The conductivity for < fj. < 
T has no analytic form as the integral in Eq. (|36[) 
involves jT^y de/\e\ 2 , which diverges as fi < T. This van- 
ishing resistivity for fi < T for impurities preserving the 
lattice symmetry is due to the fact that the scattering 
rate goes to zero faster than the chemical potential at 
the node. Thus the contribution to scattering near the 
node is dominated by other sources of scattering as com- 
pared to exchange scattering of impurities that preserve 
the lattice symmetry. 

For /3fi < 27r, fi/2n < T < fi, and T > T k we use 
Eq. ([3"0|) and Ea.(|34p and use the same approximation 
scheme as above we obtain the resistivity as 
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FIG. 4: The Kondo temperature as a function of chemical 
potential for various values of J. Energy cutoff A = 3.5eV 



and t — 2.7eV gives |J C | 
preserving case. 



9V2t 4 
4A 3 



3.94eV for the symmetry 



In above we have used the leading order correction as 
(1//3) 3 since its prefactor is fi/z which diverges as we 
take z — > 0. Higher order expansion in fi/z shows it as 
a sum of an infinite series in power of (~/i/z) n with n 
being some integer. Thus the infinite sum gives a factor 
of-1. 

In the opposite limit where /3/x ^> 27r we use Eq. (|31"j) 



Thus even for larger chemical potential the resistivity 
does not scale solely with Kondo temperature. The be- 
havior for T < (i goes like n 2 /(T 3 - Tf). 

For large chemical potential or /3[i 3> 1 we approximate 
-df/de ~ 5(e-/i). Combining with Eq.® and Eq.fgS]) 
we get 



Tk 
T 



(38) 



Thus at large chemical potential the Kondo contribution 
to resistivity is similar in the form as the magnetic im- 
purities in conventional metal. 



4J 



ttA 3 

7f 



+ -A/i 2 + (1 



-^■T k ~ — exp c 2 — 
2tt \ VM 



where C2 = 



3ir A 



3tt ;p V 2tt 



14.69. 




(35) 



1.74. 



^(3^-8) 

Thus for both cases we obtain results similar to mean 
field results obtained in Ref0- 



IV. 



From 



RESISTIVITY DUE TO IMPURITY 
CHARGE FLUCTUATIONS 



v 2 \u\ 



V 2 /\» 



e d \, for large 
JA 



~ \{€ d -n){€ d +V-n)\ 

Coulomb repulsion U, it follows that to obtain \J\ > 
the impurity level Cd must be close to the Fermi surface 
/i. Since the density of state in the graphene is propor- 
tional to the energy scale away from this Fermi surface, or 



8 



p g (e) oc |e| with p g (e) denoting graphene density of state, 
the phase space for charge fluctuation is very small and 
the local moment region is large compared with the case 
of the magnetic impurities in the conventional metal^. 
However it is still likely to have impurity level close to 
fi which is not in the local moment regiorA Thus it is 
worthwhile to estimate the resistivity contribution from 
impurity charge fluctuations. 

We use mean field approach on the Anderson im- 
purity model shown in Eq.@ and rewrite Hjj — ¥ 
^ a Un- s d\d s — Un^n^ with n s = (dld s ) determined 
self consistently, to obtain the impurity Green's func- 
tion. From the imaginary part of this Green's function 
we obtain temperature dependence of the linear response 
resistivity by assuming Boltzmann transport. Under this 
mean field approach we obtain the retarded impurity 
Green's function as£ 



n R 

^dd,. 



M = 



1 



e d -Un- 



yR 



(w) + i0+ 



(39) 



The self energy part s (w) is given by 



N 



- — E 



\e q \ 2 (oj + v) 



(40) 



q , n 



N ^-rf (uj + /.i) 2 — Vp\q\ 2 + i0 + sign(uj + p) 



mainly interested in the leading order temperature de- 
pendence of the resistivity contributed by the charge fluc- 
tuation in the Anderson impurity model. Thus we use 
the same approximation — J g(e)df/de ~ g(e)/2T 
in Eq.(J2U) to extract the leading order in temperature 
dependence. The resistivity obtained for ~ p < T is 



Psb(T) 



r{n,T) = 3T 2 a\n 



9 ni V 4 h 1 
t 2 e 2 r(/i, T) 
A 2 - I 2 



rp2 



} (3a In 



A 2 



(42) 
+ 4a 



-6) + T 2 (a ((8 + 9tt 2 ) a - 12) + 9) + 27(i d + p) 2 



Here a = ~ 2 t ^ and e d = + Un- S . Analytic result of 
resistivity for T < p can also be obtained but the ex- 
pression are cumbersome and we defer a numerical anal- 
ysis to section V. From Eq. fill) we find the non-magnetic 
region® by demanding n s = n_ s when ed — p- Within 
this charge fluctuation regime (0 < n s <C 0.5) we study 
the temperature variation of resistivity at a given p and 
Coulomb repulsion U. 



B. Czv symmetry preserving impurities 



For impurities preserving the honeycomb lattice sym- 
metry the self energy s (oj) is given by 



In above we have used |0<f, n | = |6q| 2 - |@ 9 | 2 = 1/2 for 
symmetry breaking case and |6 g | 2 = 9|g| 2 ao/8 for sym- 
metry preserving case. We take the principal part of 
Y,j ds (uj) between (—A, A) with A being the linear spec- 
trum cutoff. In the non-magnetic mixed valence regime, 
of which we are interested in, < n s = n- s < 1/2. The 
impurity occupation n s is given by 



" d^^[G B dd »] 
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7T 



(41) 



By using Eq. (|39|) and Eq.(l40f we find the relation be- 
tween td and n s by solving self-consistent conditions nu- 
merically. 



A. Csv symmetry breaking impurities 



For impurities breaking the symmetry the self energy 
H ( 

dd.s V 



^dd s( w ) obtained from Eq. (|4T)|) is given by 



J dd,. 



[UJ 



2V 2 [, , n f\(uj + p) 2 -A 2 | 
(w + fi) In ' 



9?rf 2 
+ iir\uj + p 



(w + pf 



Since T M (w) = £ s V 2 \e k \ 2 Gf ds {uj) we use Eq.® and 
Eq. ([M)l to obtain the impurity conductivity, denoted as 
v c sb( t )- The resistivity p c SB {T) = l/a c SB {T). We are 



dd.s 



2V 



2 r 



9nt 4 



(u + pfP 



xdx 
x-l 



+ ITt\uJ + p\ 



(43) 



We again use Eq.(fT8"|) and Eq. (l2"4")l to obtain the impu- 
rity conductivity, denoted as <Jg P (T). The resistivity 
Pg P (T) = l/cfgp(T). For temperature dependence we 

use the approximation — J g(e)df/de ~ 9( e )/2T in 



Eq. (f24j) to extract the leading order. To perform this 
computation we need to find the uj dependence in the 
principal integral of Eq. (|43[) . This is done by fitting 
numerically the principal value of the integral for large 
A/(uj+p). This is because the relevant integration region 
for uj in the expression of conductivity is uj C (— T, T) 
which makes A |cj + p\ in our discussion. From the 
numerical fit with 10 < fc^p^p < 10 2 (chosen for experi- 
mentally accessible range) we have 



xdx 



io x - 1 

The conductivity obtained is 



2.589 + 1.022 
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For < p < T the resistivity p c SP (T) = l/a% p (T) -> 
similar to the case for spin fluctuation in Eq. (|5o]) . For 
T < p we have 



Psp(T) 



27na h 3nTa(T 

— 7, o n I 
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t 2 r c SP (p,T) 
T((a 4 + 7r 2 )a 2 T 12 /t 4 + ...) 

3a 2 |(e e2 + M )(M 2 -T 2 ) 5 M^) 



(44) 



Here a = 2.589 is the parameter from principal integral of 
the dot self energy. We find the non-magnetic region from 
Eq. (|4ip and study the temperature variation of resistivity 
within this charge fluctuation regime (0 < n s -C 0.5). 



C. Range of validity for mean field result 
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FIG. 5: Comparison of Bethe Ansatz with mean field d level 
occupation vs impurity level for conventional metal. Fermi 
energy is set at ed = 0. Lines denoted BA means Bethe Ansatz 
results and lines denoted MF means mean field results. The 
two match better in the region where rids < 0.1. This upper 
bound increases with decreasing 17/ A. 



Before we proceed to compare the temperature de- 
pendences of resistivity due to charge fluctuations, spin 
fluctuations, and the influence of impurities position, we 
pause here to discuss the regimes where mean field results 
are valid in this Anderson impurity model. The order pa- 
rameter of this unrestricted Hatree Fock is the d level oc- 
cupation for a given spin rid,s- To make comparison with 
exact d level occupation done by Bethe Ansata 2 ^ we need 
to go back to the case for conventional metal where the 
mean field results were done by P. W. Anderson^ 3 -. The 
d level occupation for a given spin s at zero temperature 
is 



rid, 



1 

= — cot 

7T 



e<i + Un d - S - e F 



(45) 



Here A = -K~V 2 p(ep) with p(ejr) = 1/2tv the density of 
state for conventional metal. We solve Eq. (|4"5|) in the 
non-magnetic region where rid.s = n>d,-s and compare 
the answers with exact results obtained by Bethe Ansatz. 
The comparison for d level occupation for a given spin 
v.s. impurity level, with ep = 0, is shown in Figj5] From 
this figure we can see that the mean field results devi- 
ate from exact ones in a range rid, s ~ 0.06 — 0.1 for the 
range of U ~ 2 — 8A, indicating that mean field is a good 
approximation when the impurity level ea is higher than 
the Fermi energy ep or, in the other words, the impurity 
is nearby the empty orbital region. For larger Coulomb 
repulsion U/A the minimum value of of overlapping 
region is closer to ep. Since for two dimensional system 
the mean field results are marginal, we expect the mean 
field result work for €d — ep = as the s wave scat- 
tering in the conventional metal considered abov o 23 ' 26 is 
a one dimensional problem. Since the crossover shifts 
to lower and lower values of rid as U increases, this is a 
rough criterion but establishes a basis for the mean field 
calculations. 



V. RESISTIVITY TEMPERATURE 
DEPENDENCE 

In sections III and IV we have shown the analytic re- 
sults of temperature dependence of resistivity for /3fi 3> 1 
and C 1 Here we compute numerically the temper- 
ature dependence of resistivity due to spin fluctuations, 
Psb(T) and p s g P (T) for impurities breaking/preserving 
honeycomb lattice symmetry, and the temperature de- 
pendence of resistivity due to charge fluctuations, p c SB (T) 
and p% P {T). We use the full form of Rsb{w) and Rsp(ui) 
and extract the results for T > Tk with Tk obtained nu- 
merically the same way as we obtain the Kondo tem- 
perature in FigJ5] and FigfJl We compare the resistivity 
for different symmetry with the same sets of parameters. 
The resistivity due to symmetry preserving impurities is 
much smaller than that of symmetry breaking case due 
to the factor of {p/tf (see Eq.® and Eq.®). We 
examine the resistivity due to impurities spin and charge 
fluctuations in the symmetry breaking case and make 
comparison with the experimental results^ in the next 
section. 



A. Comparison of resistivity due to spin and 
charge fluctuations with different symmetry 

We use t = 2.7eV, A = 3.5eF, V = leV, and U = 4eV 
in all of the numerical results within this section. We 
choose different impurity level ed to explore the resistiv- 
ity due to spin and charge fluctuations. The resistiv- 
ity v.s. temperature is evaluated numerically between 

Tc (io^ 4 ,io- 1 ) e y. 

Let us first study the local moment region. We choose 
6d = —leV to ensure the d level occupation rid, s — 0.5. 
The chemical potential p is chosen between 10 -4 eV to 



10 




T(eV) T(eV) 



FIG. 6: Perturbation results for resistivity v.s. temperature 
for symmetry breaking case with e d = — 1 and various [i. 
The exchange coupling strength \J\ is larger than 0.2432 (for 
ji = 0.3eV) and less than 0.288461/ (for fi = 10 _4 el/) in the 
range chosen. In all cases the Kondo temperature is less than 
W~ 12 eV so the temperature range chosen is well above the 
Kondo scale. 
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FIG. 7: Perturbation results for resistivity v.s. temperature 
for symmetry preserving case with = — 1 and various /x. 
The exchange coupling strength |J| is larger than 0.2889 (for 
ji = 0.3eV) and less than 0.338061/ (for fi = 10~ 2 el/) in the 
range chosen. In all cases the Kondo temperature is less than 
W~ 12 eV so the temperature range chosen is well above the 
Kondo scale. 



0.3eV. From |J| 



v 2 \u\ 



; rr this choice of param- 

eters renders the exchange coupling strength l.lieV < 
\J\ < 1.33eV. For both cases these exchange coupling 
strengths are less than the critical value |J C | and the 
Kondo temperature obtained for both cases are ex- 
tremely small (Tfc < 10~ 12 eV). With this choice of pa- 
rameters Tk <C fJ- and therefore the analytic expression 
for Tfc corresponds to Eq. (I23|) for symmetry breaking case 
and Eq. (|35J) for symmetry preserving case. The resistiv- 
ity v.s. temperature are plotted in Fig|5] and FigJTJ 

In FigJS] we see the Tj, <§; T tails of the logarithmic 
upturns occurring when T ~ Tfc. The resistivity goes 



FIG. 8: Mean field results for resistivity v.s. temperature for 
symmetry breaking case with = and various fi. The d 
level occupation nd, s is greater than 0.024 (when \x = lQ~ 2 eV) 
and lesser than 0.093 (when fi = 0.3eV) in the chosen range. 




FIG. 9: Mean field results for conductivity at T = 10~ i eV v.s. 
chemical potential p (in unit of eV). Dots are the numerical 
data of conductivity at various chemical potential and the line 
is the fitted quadratic curve with a ~ 0.32 + 2.08/x + 51.86/x 2 . 
This quadratic behavior is also seen in Fig.3d of RefQ. 



down as chemical potential increases. This tendency 
is quite different from the case of symmetry preserving 
ones, shown in FigJTJ The dependence of resistivity on 
chemical potential fj, for symmetry preserving case shows 
p s gp(/J,) oc [i 2 by comparing the resistivity at T = lO~ 4 eV 
in Fig(7] At temperature higher than the chemical poten- 
tial the resistivity goes down with increasing temperature 
faster than the logarithmic tail for all cases in FigJTJ This 
is due to the divergence in conductivity when fi ~ T. The 
order of magnitude of resistivity at the same temperature 
for symmetry preserving case is much smaller than the re- 
sistivity for symmetry breaking case. Thus we can safely 
ignore the contributions from symmetry preserving type 
of impurities when considering the resistivity due to spin 
fluctuations. 

For the case of charge fluctuations we choose Cd = OeV 
to ensure the d level occupation nd :S < 0.1. The chem- 
ical potential \i is chosen between 10~ 4 el/ to 0.3el/. 
We compute the resistivity v.s. temperature for T C 
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FIG. 10: Mean field results for resistivity v.s. temperature for 
symmetry preserving case with td = and various fi. The d 
level occupation nd, s is greater than 0.040 (when fi — O.leV) 
and lesser than 0.093 (when n — 0.3eV) in the chosen range. 
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TABLE I: Relation between fi and td obtained by fitting with 
experimental Kondo temperature. The left two columns show 
experimental Kondo temperature at a given gate voltage. We 
compute corresponding chemical potential in the third column 
by using V g = 5.3+515.387/i 2 . The impurity level td, shown in 
the last column, is obtained by evaluating the corresponding 
exchange coupling stregth J. 



VI. COMPARISON WITH EXPERIMENTAL 
DATA 



(10~ 4 , 10 _1 )eV numerically from the mean field results. 
The resistivity v.s. temperature are plotted in Figj8]and 
Fig llOl for symmetry breaking and symmetry preserving 
cases. 

FigH shows p c SB {T) oc ln(T) for T > l(T 2 eV and 
tends to a fiat region for small temperature, which is 
very similar to the screening result of Kondo effect at 
T < Tfe. Conductivity at T — 10~ 4 shows quadratic 
chemical potential dependence, shown in FigJHl consis- 
tent with the gate voltage dependence on conductivity 
seen in the experiment 1 -. The experimental fit in Rcf.l 
for Kondo scale, however, is about one order of magni- 
tude smaller compared with the energy scale obtained in 
logarithmic temperature range in Fig|8j The tempera- 
ture dependence of resistivity in this charge fluctuation 
regime is similar to that of Kondo model in this case 
but the physics is not related to spin but charge fluctu- 
ation. To facilitate comparing our results with experi- 
mental ones in Refd we refer to the energy scale as a 
Kondo-like temperature in the following discussion. 

FigHU] also shows p c SP (T) oc ln(T) for T < [i with 
shorter range of temperature and similarly tends to a flat 
region for small temperature. The resistivity increases 
with increasing chemical potential in Fig llOl similar to 
the case of spin fluctuation. The resistivity for symme- 
try preserving case is much smaller than that for symme- 
try breaking impurities and thus ignore the contribution 
from symmetry preserving impurities. 

In summary when both types of impurities are present 
the resistivity due to impurities preserving lattice sym- 
metry is much smaller than that from impurities breaking 
the symmetry. Thus we focus our discussions on symme- 
try breaking cases for spin and charge fluctuations in the 
next section. 



Here we make comparisons with the experimental data 
given in RefyJ. We start with the perturbative results of 
Kondo model in the case of impurities breaking the hon- 
eycomb symmetry. To have large Kondo temperature 
(30K < Tfe < 90K in the experiment) the exchange cou- 
pling | J\ must be very close to its critical value | J c |. As 
perturbation breaks down when T ~ Tk, we can only an- 
alyze the gate voltage dependence of Kondo temperature 
shown in Fig. 4 of Ref|lj. The strategy is the following: We 
find the impurity level td at a given chemical potential // 
by using the experimental Kondo temperature Tk as the 
Kondo temperature obtained by the pole of resistivity, or 
1 + JRsb(0) = where J ~ V 2 /(t d - fi)(e d + U-/i). 

In the experiment the Kondo temperature is obtained 
as a function of gate voltage. We assume the gate voltage 
V g is connected with chemical potential /i via capacitive 
effect, i.e. Q/e = c g V g /e = 8c 5 /i 2 /(27 v / 37rt 2 ag) with 
Q representing the electric charges, e = 1.6 x 10~ 19 C, 
and c g = 1.15 x l0~ 8 F/cm 2 as the capacitance of the 
graphene. In the experiment of RefjT] V g = 5.3V is re- 
garded as the position of the Dirac node. Thus we take 
V g = 5.3 + 8e/ {27 V3nt 2 alc g )n 2 = 5.3 + 515. 387^ 2 by 
fixing V g = 5.3V at fx = OeV. Using the experimen- 
tal Kondo temperature at a given chemical potential we 
compute the corresponding exchange coupling strength 
J and thus determine the relationship between /i and 
impurity level t c i- The results are shown in TableJT] 

From TableU we find td oc [i. The obtained impurity 
level td changes linearly with the chemical potential \x as 
shown in lower left of Fig llll One of the main conclusions 
of this work is that the observed upturn in resistivity^ can 
be understood in terms of an Anderson impurity model 
only if the impurity level varies with the applied volt- 
age. By using the linear fit in this figure we obtain the 
Kondo temperature as a function of chemical potential 
shown in the top of FigQT] Between fi — 0.02eV to 
0.3el^ the Kondo temperature grows monotonically from 
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FIG. 11: Top:Kondo temperature as a function of chemical 
potential fi. The relationship is obtained by \J\ ~ V' 2 /]/! — 
ed | and ed oc fi which we find by fitting experimental data 
shown in the lower left figure. Lower left: Impurity level 
ed as a function of chemical potential. Dots are obtained 
by the data in Fig. 4 of RefQ] shown in Table|I] Red dashed 
line is the linear fitting function which gives e c i = —0.2254 + 
0.8951/i. Lower right: Impurity occupation of a given spin 
function of chemical potential evaluated by mean 
field approach. 




FIG. 12: Kondo temperature Tk as a function of chemical 
potential p. The relationship is obtained by |J| ~ V 2 /\p — ea\ 
and e d = -0.2254 + 0.8951/i for < p < 0.2leV and e d = 
-0.029456 (corresponds to V g = 30V or p = 0.21896^) for 
0.21eV <p< O.SeV 



14K to 90K. The decrease of Tk with increasing p for 
OeV < p < 0.02eV^ may indicate the failure of linearity 
between p and for the onset of nonzero chemical po- 
tential or the failure of the Kondo physics near the node. 
The chemical potential dependence shown in top figure 
of Fig[TT] is roughly consistent with Fig. 4 in RefQ] in the 
intermediate gate voltage. 



For gate voltage V g larger than 301^ in RefQ, the exper- 
imental Tfc begins to decrease with increasing gate volt- 
age. This can be accounted for qualitatively, as shown 
in Fig |121 by assuming that the energy of the impurity 
level no longer changes with the external gate voltage for 
V g > 30V due to sufficient charge screening. For small 
gate voltage (chemical potential close to the node) the 
experimental Tk increases monotonically with increasing 
gate voltage. In this region neither constant impurity 
level nor e<j oc p gives the corresponding experimental 
dependence on V g based on our perturbative Kondo re- 
sults. 

We also compute the impurity occupation as a func- 
tion of p by using mean field as shown by Eq. (|41[) in sym- 
metry breaking case. The obtained impurity occupation 
for a given spin increase from 0.057 to 0.082 monotoni- 
cally between p C (0, 0.3)eV. Given that the validity of 
the mean field is limited to small values of the impurity 
level occupations (see section IVC ), we expect devia- 
tions away from the mean field. Thus the system is not 
likely to stay in the local moment region near the node, 
suggesting a cross over of impurity occupation from lo- 
cal moment to empty orbital regime based on numerical 
renormalization group results in RefQJ. Thus Kondo ef- 
fect alone would not be able to explain the logarithmic 
temperature dependence seen in RefG|. 

Let us now investigate whether charge fluctuations can 
give rise to temperature dependence of resistivity seen in 
the experiment. We take ed — —0.2254 + 0.8951/x and 
evaluate the resistivity v.s. temperature from mean field 
results of impurity Green's function for symmetry break- 
ing case. For chemical potential close to the node, we 
get reasonable temperature scale (the logarithmic behav- 
ior shows up at T ~ 10~ 3 eV) from charge fluctuations 
as shown in the top figure of Fig 02] We also have p? 
being proportional to conductivity at zero temperature, 
as seen in the lower left of Fig JT3] which was observed in 
the RefQ. Rescaling p(T) by p(0) and T by Kondo-like 
temperature Tj: obtained by the temperature at which 
the resistivity begin to show logarithmic dependence in 
T, we obtain the universal curve shown in the lower right 
of Fig |13l In this range of chemical potential the impu- 
rity occupation nd, s — 0.057. It shows that even for the 
chemical potential fj, close to the node the one parameter 
scaling is still possible in this charge fluctuation scheme, 
while it is shown analytically in Eq. (|25p the one param- 
eter scaling is unlikely for \x ~ in the Kondo case. 

For chemical potential away from the Dirac node we 
plot the resistivity v.s. temperature for fi = 0.1, 0.2 and 
0.3el^ in FigfT4l The overall feature is very similar to the 
Kondo results: near zero temperature the resistivity de- 
creases with T 2 while at large temperature p(T) oc m(T). 
At zero temperature the conductivity is proportional to 
p? as in the case shown in FigJSJ However the logarith- 
mic behavior shows up at T ~ 10 _2 el^ which is about one 
order of magnitude larger than the experimental results 
in RefQ"]. Thus the charge fluctuation cannot explain the 
experiment for large chemical potential. 
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FIG. 13: Top: Resistivity v.s. temperature for charge fluc- 
tuation case with = —0.2254 + 0.8951/i and chemical po- 
tential between 10~ 4 eV and 10 _2 eV. Lower left: Quadratic 
dependence on fi for conductivity at zero temperature. Red 
line is the quadratic fitting which gives ag B (fj,,T — 0) = 
0.0000693 + 0.1675^ + 90. 3726/? and blue dots are conductiv- 
ity at various chemical potential. Lower right: Universality 
curve after rescaling resistivity and Kondo-like temperature 
scale T£ obtained by the temperature scale when logarithmic 
behavior is shown. 1.17)? ~ Tk by comparing this universality 
curve with Fig.2b in RefU. 
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FIG. 14: Resistivity due to charge fluctuations v.s. temper- 
ature plot. For T < 10~ 2 eV the resistivity decreases slowly 
with temperature. For 10 _2 eV < T the logarithmic depen- 
dence on temperature begins to appear. 
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MeV) 

FIG. 15: Kondo temperature Tk as a function of chemical 
potential \i. Blue dots are the Kondo-like temperature, ob- 
tained by comparing universal curve with NRG results, for 
charge fluctuations case. Purple dots are the Kondo tem- 
perature for spin fluctuation case. Between fi — 10 _3 eV to 
fi — 0.02eV both cases give Kondo temperature higher than 
their neighboring chemical potentials. The "mixed" (brown 
dots) case takes Kondo scale obtained by charge fluctuations 
with fi < 10~ 3 eV and spin fluctuations with n > 0.02eV. By 
combing these two the Kondo scale obtained grows monoton- 
ically with chemical potential. 



observe that charge fluctuations give large Kondo scale 
with increasing chemical potential and a good agreement 
with experimental results is obtained only if \i < 10~ 3 et^. 
Away from the node the Kondo scale obtained by charge 
fluctuations grows much faster than the that of the spin 
fluctuations. The Kondo scale obtained from spin fluc- 
tuation, on the other hand, gives large Kondo scale for 
H < 10 eV and reaches its minimum when \i ~ 0.02eV. 
The combined picture of the two cases as shown in FigfT5| 
by assuming charge fluctuation for /i ~ and spin fluctu- 
ation for large (A, can give the overall consistent picture 
as seen in the experiment for gate voltage less than 30V. 
For 10~ 3 eV < fi < 0.06eV it shows the cross over from 
charge fluctuations to spin fluctuations, which is not ac- 
counted for in our simple mean field in Anderson model 
nor perturbation in Kondo model. For gate voltage larger 
than 30V a non-monotonic dependence of T k on fi seen in 
RefQ. Given our analysis we speculate that the screening 
due to the finite density of carriers could modify the de- 
pendence of the energy of the impurity level on the gate 
voltage. A weaker dependence at large gate voltage will 
lead to a decreasing Kondo temperature. 



VII. CONCLUSION 



By comparing the universal curve obtained by numeri- 
cal renormalization group^i shown in Fig. 2b in RefQ and 
the one we have for charge fluctuations in lower right of 
FigCTJwe get 1.1T£ - T k . The "Kondo temperature" 
for both charge and spin fluctuations as a function of 
chemical potential is shown in FigfTSl From FigJT5] we 



We use Anderson impurity model to describe the di- 
lute impurities behavior in graphene. The goal is to 
test whether the recent experiment on the resistivity of 
graphene with vacancies induced by ion irradiation in 
ultra-high vacuum! can be solely explained by the single 
impurity Kondo effect (spin fluctuations). To study this 
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local moment regime we use Schrieffer Wolf transforma- 
tion to freeze the charge degree of freedom and obtain 
the Kondo Hamiltonian. In the case of dilute impuri- 
ties we may ignore the RKKY interactions and treat the 
problem as single impurity Kondo model. 

We have computed this Kondo contribution to DC re- 
sistivity by perturbation in T-matrix formulation. Ana- 
lytic expressions are obtained for 0fi <C 2ir and 0\i 2tt 
by taking asymptotic form of digamma function in the in- 
tegrand. The Kondo temperature dependence on chem- 
ical potential and exchange coupling are obtained. De- 
pending on the location of the impurities the Kondo con- 
tribution to resistivity is very different. For the type 
of magnetic impurities which break the C% v symmetry 
at low chemical potential it shows power law tempera- 
ture dependence as 1/(T — Tfc). For the magnetic im- 
purities preserving the C% v symmetry of the lattice at 
low chemical potential it shows power law dependence 
as 1/(T 3 —Tj}). At even lower chemical potential when 
Fermi surface is close to the node, both cases show ex- 
tra dependence on chemical potential as well as Kondo 
scale. Near the node a critical value of exchange coupling 
is needed for Kondo effect to be realize d 13 ! 15 . The critical 
value | J c | is larger for symmetry breaking case. For large 
chemical potential fi both cases show logarithmic depen- 
dence on temperature scaled by the Kondo temperature. 
With increasing /j, the resistivity at a given temperature 
decreases for impurities breaking the honeycomb symme- 
try while the resistivity increases for the ones preserving 
the symmetry. The resistivity obtained with same set of 
parameters show that the dominant source of resistivity 
is from the impurities which break the C^ v symmetry. 

We also have computed the effect of charge fluctuation 
for impurity occupation < 7i<2, s -C 0.5 by using mean 
field approach on the Anderson impurity model. The re- 
sistivity at a given temperature has similar dependence 
on the chemical potential as the case for spin fluctua- 
tions. Similar to the spin fluctuation case the dominant 
contribution to resistivity at the same sets of parameters 
comes from the impurities which break the honeycomb 
lattice symmetry. 

By studying the resistivity v.s. temperature and com- 
paring with experimental results in RefJT] from both spin 
and charge fluctuations in the symmetry breaking case 
we find that the Kondo effect fails to give the correct 
Kondo scale and unable to describe single parameter scal- 
ing for chemical potential nearby the node. For fi ~ 
the resistivity due to charge fluctuations give reasonable 
temperature dependence and the resistivity after rescal- 
ing also shows single parameter universal behavior. The 
same analysis yields large Kondo scale for 10~ 2 eV < fi 
in the charge fluctuation case which is roughly the same 
chemical potential at which we get non monotonic behav- 
ior of Kondo temperature in the spin fluctuation (Kondo) 
case. 

The failure of Kondo explanation nearby the node 
is consistent with the numerical RG results 14 which 
find the Kondo effect near the node is suppressed for 



r > 1/2 for systems having electronic density of state 
p(e) cx \e\ r . By combining the low chemical potential re- 
sults (0 < fi < 10 _3 ey) from charge fluctuation with the 
large 10~ 2 eV < [i results from Kondo effect obtain the 
Kondo scale consistent with the experimental results. For 
chemical potential in between these two cases the system 
should be in the mixed valence regime. For gate voltage 
higher than 30T^ a weaker dependence of the impurity 
energy on the applied gate voltage as compared to the 
dependence at smaller chemical potentials will lead to a 
decrease in the Kondo temperature. Whether this effect 
or the effect of RKKY interactions is responsible for the 
observed non-monotonic behavior on gate voltage will be 
the subject of future studies. 
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Appendix A: Derivation for symmetry breaking case 



Define y — e — /i and use 

tanh — = r »b + t— ) - V> - " 

2 iiT 2 27T A Zir 

. As tp{^ i ) h as poles on upper/lower complex plane 
we may separate Rsb(z) into two parts as Rsb(z) = 
(Ii + /2)4/(97rf 2 ) with Ii and I 2 given by 

r 1 f°° \y + iA A 2 ,,i , Jv, 



i = zi [°° a + ^ 

J-oo y y-z (2/ + /i) 2 + A 2 ^2 "2tt' 



We may write I\ — In + I\2 with 



™ J-oo y y-z{y + nY + V y 2 ' ' 2vr 
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1 .0ix 0A 
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idx 
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V>U - 1 t, 



A 
27 



idx 



-A-6 



ix - (z + fi) i(x + A) 2 2n 2 
ix 1 1 0\i 0x 

ix - (z + /x) i(x + A) ^ 2 *27r 2n 



y—z y + fi — iA 
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In the third and fourth lines of the above equation we 
replaced y by y = — jj, + ix and S — > +0. Similar compu- 
tation gives I12 as 



r A 

/l2 - Y 



+ 



dy 



y - z y + jjb + ih y + fi — iA' 



( -A) 

V 2tt> 
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j 2tt 2tx i 
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ITT 
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A iA 1 .fly pA 

2?r (/x + z) -MA W 2 l 2ir 2ir ' 



2?r J_^ y- z y + fi-iA 
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122 



dy 



V + fJ- A 2 1 ./3y 
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-(- 



■ z*y + n 
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-1 
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A . z + n 1 ,/3z 
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A 
27 

A 
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A 

ITT— — 

2ir iA 
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C 2 2 
r A+8 



z z + /i + iA 



a-5 ^ + i(fi + z))(x - A) 



dx 



iA 



1 ./3/i /3A 

07T^ ( 2 + ^ + 2^ 



z + /U 



A n 

27T Z + /i 



,,1 



The sum of 7i and I 2 is then given by 



Combining In and ii2 we get 



h = 7T- 



d _z + M ( V4 



+ i^l 

2 T * 2tt > 



z 

-A-« 



Z Z + /i 



_ 
W-oo ./-A+5' 



fi.T 



iA' 
2^(2" 



firr 
2tt 



A+5 J (x + A)(z + n - ix) 



Here Ci denotes the integration path taken from z — 
—fx — ioo to z = — jj, along imaginary axis. We may also 
write I2 = hi + I22 in different region as 



Rs B (z) = ^(h+I 2 ) 



4 rA 



97rf 2 I 7T V z — Z Z + (1 — iA 
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Rewrite z = — \i + ix in the expression of Rsb(z) along 
the C\ and C 2 paths we get 



J21 



2TT L, 



y-z (y + M ) 2 + A 2n 2 2tt ; 

y + M vii^if) 

y-z y + n + iA y + \x — iA 



z + ^(l-i§) 



2ir Jc 2 dZ z - z : + n +- iA 
A ■ " A - s 



A+8 



dx 



x + i(p + z)(x — A) 



A 



iA 



~ m YJk—{YY7)^\ + l Y + Y ) 



Rsb(z) 



1 tLti 4- §£. ) 
1 2tt ' 2ir I 



-A/zA( r dx ^L_ 

9irt 2 I in V J Q (fj, + z - ix)(x + A) 
Jo 



(H + z + ix)(x + A) J 

Z7T 







xip(± 



I A+8 



(/i + z — ia;)(a; + A) 



"27r 



2tt • 



(/i + z + ia;)(.T + A) 

-2A 2 (z + ^) 



+ ^ l 2^ 



(z + M ) 2 +A 2 



By defining F(x, ji, z) as 



Here C 2 denotes the integration path taken from z = —\i 
to z = — /U + ioo along imaginary axis. -Z22 is expressed 
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we may simplify above expression as 
Usui:-) = i d, 



j-iU? r 'dx F{x ^' z)x 

97rt 2 l7rV J x-A 



F(x, /i, z)x s 
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Appendix B: Derivation for symmetry preserving 
case 



Consider integrals of the form: 
4 r |e| 3 rfe 
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We take the integration regions into two parts by writing 
h = hi + I12 with 
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We combine results of /1 and J2 to get 
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